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Resumen 


Los modelos tradicionales para describir la extracción de agua por las plantas usan 
funciones mecanicistas O empíricas para representar el término de sumidero de 
la ecuación de transferencia de agua en el suelo, lo cual no permite representar 
explícitamente el efecto de la geometría del sistema de raíces sobre los patrones 
de extracción. Para eliminar esta limitante se desarrolla un modelo que simula la 
extracción, el cual considera las raíces como segmentos de frontera del dominio 
de solución e impone en éstos una condición tipo radiación para representar la 
resistencia al flujo del agua en la interfaz suelo-raíz, la cual es despreciada en 
los modelos reportados en la literatura que usan en las raíces una condición de 
frontera tipo Dirichlet. Se establece la condición tipo radiación, la cual contiene 
parámetros que representan las propiedades de las plantas, las características 
hidrodinámicas del suelo y la resistencia al flujo del agua en la interfaz suelo- 
raíz. Se obtiene una solución de elemento finito de la forma bidimensional de la 
ecuación de transferencia de agua sujeta en las raíces a la condición de radiación 
establecida; se aplica el modelo a un escenario de extracción, asumiendo valores 
medios característicos para las propiedades de las plantas y considerando un 
coeficiente de resistencia de la interfaz suelo-raíz variable. Se muestra que la 
resistencia en la interfaz suelo-raíz define el orden de magnitud del caudal de agua 
que extraen las raíces del suelo, por lo que esta variable no puede ser despreciada 
cuando se analiza el fenómeno de la extracción. 


Palabras clave: condición de frontera tipo radiación, coeficiente de resistencia de la 
interfaz suelo-raíz, potencial de presión de las raíces, sistema planta-suelo. 


Introducción 


Los primeros esfuerzos dirigidos al desarrollo 
de un criterio con base mecanicista para descri- 
bir la extracción de agua por las plantas se de- 
ben a Gradmann (1928) y van den Honert (1948); 
estos investigadores analizaron el flujo de 
agua del suelo a las raíces mediante relaciones 
análogas a las que describen el flujo de corriente 


eléctrica en un sistema de resistencias en serie. 
Sin embargo es hasta que Gardner (1960) pre- 
senta su modelo, cuando el estudio del proceso 
de extracción de agua del suelo por las plantas 
inicia su pleno desarrollo. Los modelos mate- 
máticos que se han desarrollado para predecir 
la extracción abordan su estudio a dos escalas 
diferentes de análisis: la microscópica y la 
macroscópica. 
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Los denominados modelos microscópicos 
describen la extracción de agua por las plantas, 
representando susraícescomotuboscilíndricos 
de longitud infinita que absorben agua uni- 
formemente del suelo. Bajo este supuesto, se 
han obtenido soluciones analíticas para calcu- 
lar la extracción, resolviendo la ecuación de 
transferencia de agua en suelo (Richards, 1931) 
en coordenadas cilíndricas y despreciando el 
efecto de la gravedad (Gardner, 1960). En esta 
escala de análisis se deduce que la extracción 
de agua es directamente proporcional a la 
diferencia de presiones entre el suelo y la raíz, 
e inversamente proporcional a la resistencia 
de la interfaz suelo-raíz. Sin embargo, al 
generalizar estos modelos a todo el sistema 
de raíces, se ha encontrado que sus resultados 
difieren sensiblemente de las observaciones 
realizadas en campo y laboratorio (Passioura y 
Munns, 1985). Esto se atribuye principalmen- 
te a la representación del sistema de raíces 
como un arreglo de tubos cilíndricos. 

Los modelos macroscópicos, por su parte, 
consideran todo el sistema de raíces de las 
plantas y pueden ser clasificados en dos 
grupos. El primer grupo está formado por 
aquellos modelos que integran las propieda- 
des de las raíces sobre un volumen elemental 
representativo de suelo, considerando con 
esta hipótesis que la raíz es una componente 
del suelo. En esta aproximación se representa 
el término de sumidero de la ecuación 
de transferencia o ecuación de Richards 
mediante funciones de extracción empíricas 
O .mecanicistas. Las funciones empíricas 
relacionan la extracción con la profundidad de 
las raíces, el contenido de humedad del suelo 
y la transpiración de las plantas, tratando de 
reproducir los datos medidos en campo sin 
tomar en cuenta la física del proceso (Feddes 
et al., 1976 y 1988; Vrugt et al., 2001; Dardanelli 
et al., 2004). Por su parte, las funciones 
de extracción mecanicistas presentan una 
estructura matemática basada en la ley 
de Ohm, y requieren de un conocimiento 
detallado de la distribución de la densidad de 
las raíces en el suelo, y de la estimación de las 


magnitudes de resistencia al flujo del agua del 
suelo y de las raíces (Nimah y Hanks, 1973; 
Feddes et al., 1974; Radcliffe et al., 1986). La 
mayor limitante de este enfoque macroscópico 
es que no permite representar explícitamente 
el efecto de la geometría de las raíces en los 
patrones de extracción y en la distribución de 
humedad en el suelo. 

El segundo grupo de modelos macroscópi- 
cos es aquel que toma en cuenta explícitamente 
la distribución espacial de las raíces en el suelo. 
El modelo representativo de este grupo es el 
desarrollado por Lafolie et al. (1991), basado 
en el análisis con la ecuación de Richards de 
cortes horizontales del sistema de raíces y 
en el uso de una condición de frontera tipo 
Dirichlet de potencial de presión, especificado 
en todos aquellos puntos donde se intersecta 
una raíz con el plano de corte. En este modelo, 
la mayor limitante es el uso de la condición 
tipo Dirichlet, que desprecia la evolución de 
la resistencia al flujo del agua de la interfaz 
suelo-raíz. 

En este trabajo se aborda el estudio de la 
extracción, considerando el sistema de raíces 
como segmentos de frontera del dominio de 
solución de la ecuación de Richards, para lo 
cual se propone lo siguiente: proyectar las 
raíces sobre el plano vertical; establecer la 
forma de la condición de frontera de radiación 
que debe usarse en las raíces para considerar 
explícitamente las propiedades hidrodinámi- 
cas del suelo y las propiedades de las raíces 
(absorción y resistencia al flujo del agua), y 
desarrollar un modelo de simulación basado 
en la solución de elemento finito tipo Galerkin 
de la ecuación de Richards bidimensional para 
analizar la extracción de agua del suelo por las 
raíces de las plantas. 


Ecuaciones de base 


La combinación del principio de conservación 
de masa y la ley de Darcy generalizada a 
medios porosos parcialmente saturados per- 
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mite obtener la ecuación que describe los 
procesos de transferencia de agua en el 
suelo. Para el caso de un suelo homogéneo, 
isotrópico, sin fuentes ni sumideros, esta 
ecuación diferencial puede escribirse como: 


98 


Er (1) 


V-[K(wv (8) 


donde 8 es el contenido volumétrico de 
agua [L*L*] y es una función del potencial 
de presión del agua en el suelo y, conocida 
como curva característica de humedad o 
curva de retención 0(y); t es el tiempo [TI]; 
V es el operador gradiente [L“]; K es la 
conductividad hidráulica del suelo [LT*], que 
en un suelo parcialmente saturado es función 
del contenido volumétrico de agua K(0) o del 
potencial de presión del agua en el suelo K(y); 
H es el potencial hidráulico [L], que es igual a 
la suma del potencial de presión de agua en el 
suelo y y el potencial gravitacional, asimilado 
generalmente a la coordenada z, orientada 
positivamente hacia abajo; H = y — z. O(w) y 
K(w) son conocidas como las características 
hidrodinámicas del suelo y son funciones no 
lineales. 

Para el análisis de la extracción con la 
ecuación (1) es necesario representar las pro- 
piedades hidrodinámicas del suelo, expresan- 
do el contenido volumétrico de agua como 
función del potencial de presión del agua en 
el suelo 6(w) y la conductividad hidráulica 
como función del contenido volumétrico 
K(0). Fuentes et al. (1992) recomiendan en 
estudios experimentales y de campo la curva 
de retención propuesta por van Genuchten 
(1980), considerando la restricción de Burdine 
(1953), y la curva de conductividad hidráulica 
propuesta por Brooks y Corey (1964), debido a 
que esta combinación satisface las propiedades 
integrales de la infiltración y a la facilidad para 
identificar sus parámetros. 

La relación propuesta por van Genuchten 
(1980) es: 


—m 


o(y)=0, + (8, —0,)|1 a] (2) 
y 


d 


donde 6, es el contenido volumétrico de agua 
a saturación; 6, el contenido volumétrico 
residual; y y Un valor característico de la 
presión del agua; m y n, parámetros de forma 
empíricos relacionados a través de la restric- 
ción de Burdine (1953): m=1-2/n, con0<m 
<lyn>2. 

El modelo para la conductividad hidráulica 
propuesto Brooks y Corey (1964) es: 


(3) 


n 
9-6 

K0)=K/; Se 
E r 


donde K_ es la conductividad hidráulica a 
saturación [LT*] y n es un parámetro de forma 
adimensional. 


Condición de frontera en las raíces 


La condición de frontera que se debe utilizar 
en la frontera de las raíces se deduce de la 
siguiente manera: 

Por una parte, el flujo de agua en el suelo 
en la dirección normal al segmento de frontera 
que representa las raíces (q,,,) se puede calcular 
utilizando la ley de Darcy: 


9 (y-z) 
on 


Uns = -K (y) (4) 


donde n denota normal a la frontera. 

Por otra parte, en analogía con la ley 
de Ohm, el flujo de agua en las raíces en la 
dirección normal (q, ) está dado por: 


Qnr = Kia (| mo pe | + Es (W;,) cos O, (5) 
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donde y, es el potencial de presión del agua 
en interior de la raíz; K, (y, ), la conductividad 
de la interfaz suelo-raíz [LT"]; A, una longitud 
característica de la interfaz suelo-raíz [L], y 
cosó_ es el coseno director asociado con la 
coordenada vertical. 

La resistencia de la interfaz suelo-raíz [T] 
está dada por: 


y 


6 
Kv Ko (4) / 


Rin(Wj,) = e 


Por continuidad en un punto localizado 
sobre la superficie de la raíz, el flujo de 
agua en el suelo q,.es igual al flujo de agua 
en la raíz q,,; para cumplir esta condición se 
igualan las ecuaciones (4) y (5), y al considerar 
la expresión (6) se obtiene la condición tipo 
radiación (Carslaw y Jaeger, 1959) siguiente: 


(y e 2. Y yx ., (95) 0050, (7) 


(y, 
Ria (Vi) 


K sn 


Para aplicar la ecuación (7) es necesario 
disponer de relaciones funcionales para 
los parámetros que en ella intervienen. 
Considerando resultados de estudios que 
muestran que la resistencia en la interfaz es 
definida principalmente por la raíz (Taylor 
y Klepper, 1975; Herkelrath et al., 1977), en 
una primera aproximación se propone que 
la conductividad de la interfaz suelo-raíz se 
calcule como el promedio geométrico de la 
conductividad hidráulica del suelo K(y) y la 
conductividad de la raíz, la cual es función de 
la presión del agua en la raíz K (y): 


Kia (Y; )= JK(wK, (y,) (8) 


El comportamiento extremo de la relación 
(8), considerando las ecuaciones (6) y (7), 
permite deducir que la extracción es nula 
cuando la conductividad en la interfaz tiende 
a cero; mientras que la extracción es máxima 
cuando la conductividad de la interfaz alcan- 


za su valor máximo. En consecuencia se tiene 
que la condición tipo radiación (7) describe el 
fenómeno de la extracción y contiene como 
caso particular la condición tipo Dirichlet 
propuesta por Lafolie et al. (1991), que 
considera despreciable la resistencia al flujo 
del agua en la interfaz suelo-raíz, hipótesis 
que físicamente no es correcta. 


Modelo numérico 
Se desarrolla un modelo basado en la solución 


numérica de la forma bidimensional de la 
ecuación (1): 


98 _ 0 ow| 0 ow| 0K(y) 
ol 2 K00 le Kw a 


donde x y z son las coordenadas espaciales 
horizontal y vertical, respectivamente [L]. 

La integración en el espacio de la ecuación 
(9) se realiza con el método del elemento 
finito tipo Galerkin, considerando funciones 
de interpolación lineales y su integración 
en el tiempo con un método de diferencias 
finitas implícito; se usa el método iterativo 
de Picard para linealizar el sistema de 
ecuaciones resultante de la discretización 
y un método de gradiente conjugado 
precondicionado para resolver el sistema de 
ecuaciones algebraicas (Noor y Peters, 1987), 
el cual ha sido adaptado para utilizar un 
almacenamiento matricial libre de ceros. Los 
detalles del procedimiento de discretización 
están bien documentados en la literatura 
y pueden ser consultados, por ejemplo, en 
Zavala y Ojeda (1999), o en Zienkiewicz et al. 
(2005). El esquema numérico que se obtiene al 
incluir la condición de frontera tipo radiación 
(7) es el siguiente: 


12. pH 2 
Hao Wi =N E (10) 


siendo: 
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M!; 
eE ra? 0 


CPF y pa 
o ra E 
ra —0 
mp2 AY ! (12) 
At Af 


tesel tiempo; At, el paso de tiempo); p, el número 
de iteraciones en el intervalo de tiempo; C 
= d8/dy, la capacidad específica (Richards, 
1931); M, la matriz de masa; B, la matriz de 
rigidez; G, el vector de la fuerza de gravedad; 
Q, el vector de flujos en la frontera asociados 
con una condición tipo Neumann, F y L, la 
matriz y vector asociados con la condición de 
frontera tipo radiación. Los coeficientes de 
estas matrices y vectores son: 


M5 (13) 

B, Ez + 0,0,) (14) 
=> + (15) 

O; =24, a (16) 
ES Lan 
Ea (18) 


donde e es el número de elementos finitos; Sy, 
la delta de Kronecker; A, el área del elemento 
finito triangular; K, la conductividad hidráulica 
promedio en el elemento; u y U, factores 


geométricos definidos según las funciones de 
base u, = 2, 2,Y 0,=X,—X, los subíndices 1, j, k, 
que corresponden a las esquinas del elemento 
finito triangular, corren sobre sus tres posibles 
permutaciones secuenciales (Zienkiewicz et al., 
2005); Ro el valor de la resistencia en la interfaz 
suelo-raíz en cada nodo sujeto a la condición 
tipo radiación; L, la longitud de frontera 
asociada con suda. nodo sujeto a la condición 
tipo radiación; q, un valor de flujo de agua 
en con una condición tipo Neumann, 

a, Ya longitud de frontera asociada con cada 
ea sujeto a una condición tipo Neumann, V,, 
el valor de la presión en la raíz en cada nodo 
sujeto a la condición tipo radiación. 

El caudal de agua por unidad de longitud 
que se transfiere a través de los nodos 
localizados en las fronteras sujetas a la 
condición tipo radiación (Q,) se calcula de la 
siguiente forma: 


+1/2 
0 ñas == via), Kolas cos 0, L, q) 
Rontrar 


El esquema numérico (10-19) se codificó 
en lenguaje de programación C++ y permite 
realizar modelación de escenarios de extracción 
con diferentes condiciones iniciales y de 
frontera. 


Aplicación y resultados 


Tomando en cuenta que la distribución 
bidimensional del sistema de raíces de las 
plantas puede ser caracterizada, por ejemplo, 
mediante modelos de ruptura dieléctrica 
(Chikushi e Hirota, 1998) o patrones fractales 
(Tatsumi et al., 1989), es posible generar conestos 
métodos dominios de solución para alimentar 
el modelo de extracción desarrollado. En este 
trabajo se modela la extracción considerando 
una proyección del sistema de raíces de una 
planta de mijo (Panicum miliaceum) sobre el 
plano vertical (x-z) reportada por Tatsumi et 
al. (1989); para facilitar la visualización de 
resultados, a la imagen original se le aplicó un 
factor de escala de 10:1 (figura 1). 


Tecnología y 


Ciencias del Agua, vol. UL, núm. 1, enero-marzo de 2012 


Ciencias del Agua, vol. 11, núm. 1, enero-marzo de 2012 


Zavala et al., Modelo bidimensional para la extracción de agua por las raíces de las plantas 


Figura 1. Proyección de las raíces de una planta de mijo 
(Tatsumi et al., 1989). 


La condición inicial que se considera en 
la modelación numérica corresponde a una 
distribución de humedad uniforme cercana 
a la saturación, que puede representar, por 
ejemplo, la condición de humedad que se 
alcanza instantes posteriores a la aplicación de 
un riego: 


y(x,z)=V j, t=0 (20) 


donde y, ,=-0.10 cm. 

Las condiciones de frontera consideradas 
en los segmentos del dominio distintos de 
aquellos que representan el sistema de raíces 
son (ver figura 1): 1) en la superficie del suelo 
se impone una condición de frontera tipo 
Neumamn de flujo variable para representar 
la evaporación superficial; 2) en la fronteras 
laterales se impone una condición de flujo 
nulo; 3) en la frontera inferior se usa una 
condición de gradiente unitario. Este conjunto 
de condiciones de frontera se pueden escribir 
como: 


ke. + 


Oz máxi 


e A O<x<a y 
24 


(la) 


b<x<x 


N 
11 
o 


máx” 


xy) 002. 1 OS<SX<Xmáx; Z =Zmáx (21b) 


090 0 


x=0 
Ox le 


X= Xmáx;s 


O<z<z (21c) 


máx 


donde E, ., =0.2cm/d;x,, =200cmyz,. = 
200 cm. Alo largo de la frontera ab se impone la 
condición tipo radiación (7). 

El dominio mostrado en la figura 1 se 
discretizó generando una malla de elementos 
finitos con el programa ARGUS ONE 4.0p, 
respetando la condición de uniformidad 
(Zienkiewicz et al., 2005) y concentrando los 
nudos cerca de la raíz y la superficie, que es 
donde se generan los mayores gradientes de 
presión (figura 2). La malla consta de 25773 
nudos y 45395 elementos, con espaciamientos 
mínimos: Ax. =1.38 x 10* cm y Az, . =1.38 x 
10 cm. 

Se utiliza el suelo franco de Fuentes ef 
al. (1992), cuyos parámetros hidráulicos son 
los siguientes: O. = 0.45 cm*/cm”, O, = 0.10 
cm*/cm?, y, = -100 cm, m =0.5, n = 12 y una 
conductividad hidráulica a saturación K_= 6 
x 10* cm/s. El potencial de presión en la raíz 
y la conductividad hidráulica se asumieron 
constantes, utilizándose los siguientes valores 
característicos y, = 1.0 x 10* cm y K_ =1x 10* 
cm/s (Taylor y Klepper, 1975; Radcliffe et al. 
1986), y una longitud característica de A, =1x 
10 cm. 

El paso de tiempo inicial para la simulación 
fue de At, =0.1 s, y éste se incrementó 10% si 


ini 
M4 ” 


el número de iteraciones “p” realizadas en el 
paso de tiempo era menor de 3, hasta alcanzar 
el valor máximo de Af,.. = 300 s; mientras 
que el criterio de convergencia adoptado fue 


_ 4 p+l p+1/2 y) 3 
e= máx lua, ON o <1x10 “cm. 
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Figura 2. Malla de elementos finitos. 


En las figuras 3 y 4 se presentan las líneas 
equipotenciales a los 7 y 21 días obtenidas con 
el modelo numérico; se puede observar que 
la humedad disminuye rápidamente en las 
zonas de suelo que presentan mayor densidad 
de raíces, lo cual muestra la importancia de 
representar explícitamente la arquitectura del 
sistema de raíces en el estudio de la extracción. 
En la figura 5 se presenta la variación en el 
tiempo del gasto extraído por la raíz, teniéndo- 
se que la extracción disminuye sensiblemente 
durante los primeros tres días de simulación 
y posteriormente presenta una disminución 
gradual; este comportamiento se debe a la 
evolución de la resistencia en la interfaz suelo- 
raíz. Al inicio de la simulación, la resistencia en 


la interfaz suelo-raíz alcanza su valor mínimo 
(R, = 4.1 x 10* s) y la extracción alcanza valores 
máximos conforme el contenido de humedad 
del suelo que circundan a la raíz decrece, la 
resistencia enlainterfaz suelo-raízseincrementa 
y el gasto extraído disminuye, evidenciando 
que la resistencia en la interfaz suelo-raíz define 
la magnitud de caudal de agua extraído por las 
plantas, tal y como han mostrado numerosas 
investigaciones experimentales (Passioura y 
Munns, 1985). 

La condición de frontera tipo radiación (7), 
además de incluir la resistencia al flujo del agua 
en la interfaz suelo raíz, considera también el 
potencial de presión del agua en el interior de 
las raíces, el cual determina la apertura de los 
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Figura 3. Equipotenciales a los siete días (cm). 
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Figura 5. Variación en el tiempo del gasto de agua extraído 
por ancho unitario de raíz, calculado con la ecuación (19). 


estomas, y por lo tanto la transpiración, por 
medio de mensajes químicos. En consecuencia, 
una buena descripción y caracterización de 
los parámetros de la condición de radiación 
(presión y conductividad de las raíces) permite 
un análisis más detallado de la extracción 
respecto del que puede obtener con los modelos 
reportados en la literatura que manejan 
condiciones de frontera tipo Dirichlet o tipo 
Neumamn. 


Conclusiones 


Se ha propuesto una aproximación mecanicis- 
ta para describir el fenómeno de la extracción de 
agua del suelo por las plantas, que está basada 
en el análisis con la ecuación de Richards de 
proyecciones del sistema de raíces sobre el 
plano vertical. 

Se ha establecido la condición de frontera 
tipo radiación que debe sujetar la ecuación 
de Richards en la frontera de las raíces de las 
plantas para describir la extracción de agua 
del suelo; esta condición de frontera incluye 
parámetros hidráulicos del suelo, propiedades 
de la planta y un coeficiente de resistencia 
al flujo del agua de la interfaz suelo-raíz 
que es variable, lo cual permite describir 
el comportamiento de la extracción bajo 


diferentes condiciones de humedad en el suelo 
y necesidades hídricas de la planta. El análisis 
del comportamiento extremo de la condición 
de radiación establecida ha permitido mostrar 
que contiene la condición de frontera tipo 
Dirichlet propuesta por Lafolie et al. (1991). 

Se ha desarrollado un modelo basado en la 
solución numérica de la ecuación de Richards 
bidimensional y se ha simulado un escenario 
de extracción, cuyas condiciones asimilan, por 
ejemplo, las condiciones generadas después 
de la aplicación del riego. El escenario 
de extracción analizado muestra que la 
conductividad de la interfaz suelo-raíz define 
el orden de magnitud del caudal o gasto de 
agua que extraen las raíces de las plantas y 
evidencia la importancia de caracterizar este 
parámetro. El modelo desarrollado es una 
herramienta útil para simular escenarios 
de extracción en diferentes tipos de suelo y 
cultivos. 
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Zavala et al., Modelo bidimensional para la extracción de agua por las raíces de las plantas 


Abstract 


ZAVALA, M., SAUCEDO, H., BAUTISTA-CAPETILLO, C., € FUENTES, C. Two- 
dimensional model for water uptake by root systems. Water Technology and Sciences (in 
Spanish). Vol. IL, No. 1, January-March, 2012, pp. 33-43. 


Classical models to describe water uptake by roots use mechanistic or empirical functions to 
represent the sink term in the water transfer equation, which does not allow for analyzing the 
effect of the root system geometry on extraction patterns. A model to simulate water uptake 
was developed in order to eliminate this limitation. This model considers the roots to be one 
boundary of the domain and uses a radiation condition which allows for representing the water 
flow resistance at the soil-root interface. This resistance is neglected by the models reported in 
the literature because they use a Dirichlet boundary condition for the roots (specified pressure). 
The radiation condition is established, which contains parameters that represent the properties 
of the plants and hydraulic soil as well as the resistance to water flow at the soil-root interface. A 
finite element solution is obtained for the two-dimensional form of the water transfer equation 
subject to the radiation condition at the roots. This model is applied to the analysis of a water 
uptake case, with typical values for the plants” properties and considering the variability of 
the resistance of the soil-root interface. Since the resistance at the soil-root interface is shown 
to define the magnitude of the water uptake by roots, this variable cannot be neglected when 
analyzing the water uptake phenomenon. 


Keywords: radiation boundary condition, soil-root resistance coefficient, root water potential, 
soil-plant system. 
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